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ABSTRACT 

Three  computer  simulations  of  command  to  line-of-sight 
missile  guidance  systems  were  developed  within  an  electronic 
warfare  environment  to  test  the  feasibility  of  improved 

* 

guidance  techniques  while  maintaining  simplicity  for  generic 
application  among  the  various  missiles  of  this  general 
guidance  classification.  The  first  simulation  modeled  the 
basic  guidance  scheme;  the  second  introduced  a "lead-angle" 
concept;  and  the  third  simulation  combined  these  techniques 
for  use  depending  upon  whether  the  target  took  evasive 
action  or  employed  electronic  countermeasures.  It  was  found 
that  consideration  should  be  given  for  use  of  a "lead-angle" 

* variant  in  conjunction  with  the  basic  guidance  technique  to 

enhance  the  effective  engagement  envelope  of  these  missile 
systems  against  relatively  slower-maneuvering  targets  employing 
only  low-duty-cycle  denial  jamming. 
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I . INTRODUCTION 


Command  to  line-of-sight  guidance  for  several  types  of 
surface-to-air  missile's  was  developed  through  computer 
simulation.  Consideration  was  given  to  electronic  counter- 
measures against  the  guidance  methods  of  this  type. 

The  simulations  are  of  sufficient  complexity  to  demonstrate 
feasibility  while  maintaining  simplicity  for  generic  application 
among  the  various  missiles  of  this  general  guidance  classi- 
fication. 
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II.  OVER-VIEW  OF  GUIDANCE.  TYPES 


Guidance  systems  can  be  considered  to  fall,  generally,  into 
one  of  three  basic  types:  line  of  sight,  pursuit,  or  propor- 
tional navigation. 

The  line-of-sight  system  maintains  the  missile  on  the  vector 
from  the  ground-based  point  of  missile  control  (the  tracker) 
to  the  target  aircraft.  Variations  to  this  system  can  provide 
for  "lead-angle"  intercept  through  missile  command  computer 
generation  of  a synthetic  line  of  sight  upon  which  the  missile 
is  commanded  to  fly  for  target  intercept.  Also,  beam-riding 
missiles  fit  this  category.  However,  in  the  beam-riding  type  of 
line-of-sight  guidance  the  missile  avionics  are  more  complicated 
and  expensive  because  the  missile  normally  senses  its  position 
relative  to  the  beam  center  and  generates  its  own  guidance 
commands  internally.  This  extends  the  missile's  conceivable 
range  but  usually  necessitates  some  form  of  terminal  homing. 
Basic  command  to  line-of-sight  systems  can  utilize  relatively 
cheap  missiles  by  performing  guidance  computations  in  the 
tracker  unit  and  not  in  the  missile  itself.  This  limits  the 
missile  to  short-range  scenarios  and  requires  that  the  missile 
possess  a greater  speed  advantage  over  the  target. 

Pursuit  guidance  systems  maintain  the  missile's  velocity 
vector  toward  the  target.  Generally,  sensors  in  the  missile 
can  be  assumed  to  perform  this  function;  therefore,  the  missile 
can  be  forgotten  by  the  launcher  after  the  missile  has  locked 
onto  the  target.  This  type  of  system  is  capable  of  medium 
range  work  against  a slowly-maneuvering  target.  It  has 
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relatively  simple  processing  avionics. 

Proportional  navigation  guidance  systems  hold  a constant 
angle  between  the  missile- to-target  line  of  sight  and  the 
missile  axis,  thereby,  generating  a " constant-bearing-decr easing- 
range"  , collision  situation.  In  other  words,  the  change  in 
missile  heading  is  adjusted  proportionally  to  the  rate  of 
change  in  the  missile-to-target  line  of  sight.  This  type  of 
guidance  is  good  for  long  ranges  (and  short)  against  highly- 
maneuvering  targets.  However,  the  missile  requires  either 
seeker  and  angle-rate  sensing,  or  a command  guidance  signal 
based  upon  a collision  course. 
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III.  TYPICAL  ENGAGEMENT  SEQUENCE' 


In  order  to  provide  a "vehicle"  through  which  to  better 


understand  the  basic  aspects  of  command  to  line-of-sight 


guidance  methodology,  the  engagement  sequence  of  a short-range, 


air-defense,  missile  system  is  described.  The  Roland  system 


was  selected  because  the  general  operational  aspects  of  the 


system  are  available  at  the  unclassified  level  [1] . 


The  entry  of  one  or  more  aerial  targets  into  the  range  of 


the  search  radar  is  indicated  to  the  Roland  vehicle  commander 


by  an  audible  tone.  At  the  same  time,  a synthetic  display  of 


the  targets  appears  on  a screen  to  give  the  commander  all  the 


information  needed  to  select  the  most  threatening  target.  The 


screen  images  are  different  for  friendly  and  enemy  targets. 


Also,  the  entry  of  the  target  into  the  missile  envelope. 


utilizing  target  advanced-range  computations,  is  indicated  by 


a change  in  the  display.  With  the  search  antenna  raised  and 


the  search  radar  activated,  target  acquisition  is  possible  even 


when  the  vehicle  is  in  motion. 


There  are  three  modes  of  identification,  friend  or  foe 


(IFF)  interrogation:  automatic,  manual,  and  automatic  within 


a given  range. 


When  the  commander  has  recognized  a target  as  hostile 


and  decided  to  engage  it,  he  places  a cursor  over  the  screen 


image.  This  automatically  brings  the  turret  to  bear  and 


tracking  can  commence  in  either  the  "radar"  or  "optical"  modes. 
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In  the  "radar”  mode,  the  tracking  radar  automatically 
accepts  target  designation  from  the  search  radar,  searches 
for,  locks  onto,  and  tracks  the  target. 

In  the  "optical"  mode,  the  aimer  searches  for  the  target 
in  elevation  with  an  optical  sight.  To  aid  him  an  electronic 
instrument  displays  the  maximum  theoretical  elevation  for 
the  search.  When  the  aimer  has  acquired  the  target  in  his 
cross-hairs,  he  keeps  the  target  in  his  sight  by  manipulating 
a control  stick.  This  control  keeps  the  target  properly 
positioned  by  moving  the  turret  in  azimuth  and  swivelling  a 
mirror  in  elevation. 

As  soon  as  the  commander  confirms  that  the  target  is 
within  missile  range,  he  initiates  the  firing  sequence  in 
the  "radar"  mode,  or  authorizes  "optical"  mode  firing  through 
a command  displayed  in  the  aimer's  sight.  The  aimer,  then, 
can  initiate  the  firing  sequence. 

The  missile  is  guided  by  a command  to  line-of -sight 
technique.  This  means  that  the  target  is  tracked  optically 
or  by  radar  and  the  deviation  of  the  missile  from  this  line 
of  sight  is  determined  and  corrected  by  a guidance  command. 

The  comm-'nder  may  switch  from  "radar"  to  "optical"  and  back 
again,  as  desired,  even  after  the  missile  has  been  launched. 

Target  tracking  and  determination  of  the  missile's  deviation 
from  the  line  of  sight  are  different  for  each  mode.  In  the 
"radar"  mode,  the  guidance  radar  has  two  receiving  channels. 

One  is  used  for  target  tracking  and  the  other  is  used  to 
locate  the  missile  in  the  radar  lobe  through  reception  of  the 
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missile's  radio  frequency  beacons.  By  comparing  these  angles, 
an  error  between  the  missile  and  the  target  line  of  sight  can 
be  determined.  In  the  "optical"  mode,  a biaxially-stabilized 
mirror  is  manually  controlled  to  keep  the  target  vertically 
in  the  aimer's  sight  and  the  turret  is  rotated  to  the  azimuth 
of  the  target  line  of  sight.  An  infrared  goniometer  is  mounted 
to  provide  missile  angle  from  the  tracker  by  following  flares 
mounted  on  the  rear  of  the  missile.  Then,  a deviation  of  the 
missile  angle  from  the  target  line  of  sight  can  be  determined. 

Two  groups  of  signals  are  introduced  into  the  command 
computer:  the  velocity  of  the  line  of  sight  in  azimuth  and 
elevation,  and  the  deviation  of  the  missile  from  the  line  of 
sight  in  azimuth  and  elevation.  Based  upon  data  from  the  line- 
of-sight  movement  and  the  angular  deviations  of  the  missile, 
the  necessary  guidance  signals  are  calculated. 

The  guidance  signals  are  relayed  to  the  missile  by  a 
transmitter  with  highly  directional  characteristics.  The 
command- transmitting  antenna  is  slaved  to  the  missile  angle 
in  both  azimuth  and  elevation.  It,  therefore,  is  trained 
on  the  missile  continuously. 

The  side  forces  required  for  missile  course  corrections 
are  produced  through  deflection  of  the  exhaust  jet  of  the 
sustainer  motor  by  spoilers  at  the  rear  of  the  missile  (thrust- 
vector  control) . 

When  the  missile  reaches  the  point  of  impact  with  the 
target,  the  warhead  is  detonated  by  either  percussion,  contact 
fuse  or  the  radio-frequency,  proximity  fuse.  The  warhead 
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consists  of  a • radial-effect,  multiple- fragmentation  charge. 

Figure  III-l  presents  an  operational  schematic  of  the 
basic  Roland  missile  system  operation. 
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Figure  III-l  Roland  Missile  System  Operational  Schematic 

The  computer  simulations  contained  herein  are  generic  in 
nature  within  the  command  to  line-of-sight  guided-missile  type 
and  have  only  reasonable  estimates  of  missile  capabilities 
introduced.  This  ensures  unclassified  results.  At  the  same 
time,  the  simulations  are  of  sufficient  complexity  to 
properly  weigh  the  relative  merits  of  the  guidance  variations 
discussed. 
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IV.  COMMAND  TO  LINE-OF-SIGHT  GUIDANCE  GEOMETRY 


There  are  many  tactical,  inherently  short-ranged,  surface- 


to-air  missiles  which  use  this  guidance  technique  for  their 


operation.  Its  objective,  basically,  is  to  keep  the  error 


term  (e)  zero  so  that  the  missile  flies  directly  up  the  line 


of  sight. 


Since  only  angles  to  the  target  and  to  the  missile  are 


required  for  guidance,  attempts  to  preclude  target  range 


through  denial  jamming  prove  useless  against  this  type  of 


guidance  system.  The  only  necessary  information  is,  in  general. 


to  know  whether  the  target  is  within  the  missile  envelope. 


Deceptive  down-link  jamming,  however,  appears  to  be  feasible. 


One  possible  method  is  presented  in  Section  VI. 


The  geometry  depicted  in  Figure  IV-1  summarizes  the 


"basic"  command  to  line-of-sight  geometric  situation. 
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Figure  IV-1  "Basic"  Command  to  Line-of-Sight  Geometry 


A variant  of  this  ''basic”  command  to  line-of-sight  guidance 
approach,  which  shows  promise  for  increasing  system  effectiveness, 
is  to  incorporate  "lead-angle”  into  the  basic  guidance  model. 

This  technique  would  increase  the  effective  range  of  the  missile 
against  a non-maneuvering,  crossing  target;  however,  it  is  de- 
graded by  a highly-maneuvering  target.  Since  range  information 
is  required  to  calculate  the  impact  point,  additional  electronic 
countermeasures  could  be  employed  by  the  target  to  deny  this 
range  information  thus  negating  the  "lead-angle"  command  to  line- 
of-sight  missile  guidance  scheme. 

The  previous  "basic"  geometry  has  been  retained  but  it  has 
been  modified  to  utilize  a synthetic  line  of  sight  upon  which 
the  missile  flies.  This  synthetic  line  of  sight  is  the  direction 
of  the  impact  point  from  the  ground  tracking  unit.  Therefore, 
the  missile  is  commanded  to  fly  toward  the  impact  point  rather 
than  directly  toward  the  target.  This  method,  then,  reduces 
the  missile's  flight  time  until  impact,  thus,  improving  the 
missile  system's  effective  engagement  range.  Figure  IV-2  shows 


Figure  IV-2  "Lead- Angle"  Command  to  Line-of-Sight  Geometry 


The  guidance  models  used  in  the  computer  simulations  utilize 
on-off  control.  This  type  of  control  also  is  known  as  relay 
control  or  "bang-bang"  control  because  the  system's  actuator  is 
not  moved  in  proportion  to  the  error  signal  but,  rather,  is 
positioned  either  totally  in  the  positive  direction  or  totally 
in  the  negative  direction.  This  causes  the  missile's  thrust- 
vector  control  to  apply  missile-maneuvering  thrust  fully  in  one 
direction  or  the  other.  At  first  glance  this  type  of  system 
might  appear  to  be  very  crude;  however,  it  does  have  major  advan- 
tages. First,  the  system  does  not  require  complex  proportional 
amplification;  the  unit  is  smaller  than  its  proportional  control 
counterpart;  it  is,  generally,  the  time— optimal  control  system; 


it  weighs  less  than  other  systems;  and,  finally,  it  is  relatively 


inexpensive  [2] . 
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Figure  V-l  Missile  Cross-Range  Error  with  Zero  Initial 
Values  for  Cross-Range  Error  and  Error  Rate 
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Based  upon  this  type  of  missile  control  system,  initial 
computer  simulations,  not  presented  here  in  detail,  were  designed 
They  provided  an  elementary  computer  program  structure  which  was 
modified  incrementally  to  generate  improved  computer  models.  In 
these  first  programs,  the  error  values  were  introduced  directly 
into  the  missile's  thrust-vector  control  without  error  rate. 

This  produced  increased  missile  cross-range  error  with  flight 
time  as  shown  in  Pigure  V-l.  This  effect  has  been  eliminated 
in  the  subsequent  computer  guidance  models  by  proper  acceleration 
limiting  techniques.  Additionally,  if  initial  erroneous  values 
of  cross-range  error  and  cross-range  error  rate  are  introduced 


Cross-Range  .4I)| 
Error  Rate 
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Figure  V-2  Phase  Plane  When  Control  Function  By-Passed 

With  Initial  Values  Introduced  for  Cross-Range 
Error  and  Error  Rate 
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missile  moves  about  the  line  of  sight  with  a cross-range  error 
varying  as  shown  in  Figure  V-5  - The  cross-range  error  does 
increase  by  a small  amount  with  time  as  shown  in  Figure  V-l. 

The  models  discussed  later  are  somewhat  unrealistic  because 
they  assume  that  the  required  force  obtained  from  the  missile's 
thrust-vector  control  can  be  applied  instantaneously  at  the 
proper  time  and  reversed  in  a similar  fashion  to  decelerate  the 
missile.  This,  however,  does  provide  optimum  missile  performance 
for  comparison  with  any  real  system. 

In  all  the  simulations  the  missile  is  acceleration-limited 
with  instantaneous  application  of  constant  side-force  assumed. 

One  could  allow  for  a time  delay  in  switching  the  missile 1 s 
thrust-vector  control  from  the  positive  to  the  negative  direction. 
However,  chis  would  raise  the  order  of  the  analysis.  The  computer 
program,  though,  can  be  modified  to  handle  systems  to  the 
thirtieth  order. 

Optimum  acceleration  switching  requires  that  a switching- 
boundary be  established  that  zeros  the  error  in  the  minimum  time. 
That  boundary  was  determined  to  be  parabolic  as  drawn  in  the 
phase-plane  representation  of  Figure  V-6  and  as  supported 
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Figure  V-6  Optimum  Parabolic  Switching  Boundary 
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by  tha  calculations  which  follow  later.  When  the  acceleration 
or  deceleration  trajectories  attempt  to  cross  this  optimum 
curve  (switching  boundary)  the  system  switches  to  that  optimum 
path  and  the  error  is  driven  to  zero  along  it. 

In  order  to  verify  this  result,  a phase  plane  was  produced 
by  introducing,  as  initial  conditions,  large  cross-range  error 
and  cross-range  error  rates  into  a computer-simulated  guidance 
model  containing  acceleration- limiting  methods  within  the  tracker 
control  function.  The  resulting  phase  plane,  Figure  V-7,  sub- 
stantiates that  the  tracker  control  function  switches  properly 
to  the  theoretically-optimum,  parabolic  curve  for  error  reduction 
expected  by  calculation.  This  curve  is  represented  in  Figure 
V-6.  Also,  the  target  was  intercepted  by  the  missile  as  sub- 
stantiated in  Figures  V-8  and  V-9.  Figure  V-10  shows  the 


Cross-Range 
Error  Rate 
(m/sec) 


Cross-Range  Error  (m) 


Figure  V-7  Phase  Plane  When  Control  Function  Utilized 
With  Initial  Errors  Introduced  for  Cross- 
Range  Error  and  Error  Rate 
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missile's  cross-range  error  with  time  of  flight.  This  error  is 
not  readily  apparent  in  Figure  V-9  because  of  the  relatively 
large  distance  scale  utilized  Li  that  figure.  As  can  be  seen 
this  cross-range  error  is  reduced  quickly. 

The  theoretical  calculations  supporting  the  parabolic 
nature  of  the  switching  boundary  represented  by  Figure  V-6 


follow: 


m = missile  mass 

F = lateral  force  applied  to  missile 

CRE  = cross-range  error 

CRE  = cross-range  error  rate 
• • 

CRE  = cross-range  error  acceleration 

G = magnitude  of  cross-range  missile  acceleration 
U = cross-range  missile  acceleration 
C = constant 

Since  the  mass  of  the  missile  and  the  magnitude  of 
the  output  force  from  the  missile  thrust-vector 
control  are  assumed  constant,  the  magnitude  of  the 
acceleration  is  constant,  thus: 


m CRE  = F 


CRE | = |F|/m  = G 


- CRE  I = G 


J*  | d CRE  | = Jg  dt  ->|  CRE  | = Gt  + c 

Let  the  initial  conditions  for  time  and  cross-range 
error  be  zero.  Therefore,  the  constant,  "C", 
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is  zero  also, 

CRE  = Gt  , CRE  - 0 
t = CRE/G  -*■  t2  = CRE2/G2 

| CRE  = Gt  / d (CRE)  = /Gt  dt 

CRE  = Gt2/2  + j CRE ( 0 ) | , t2  = CRE2/G2 
CRE  = CRE2/2G  + | CRE (0) J , CRE  - 0 
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Since  the  direction  of  the  missile's  thrust-vector  *: 

control  force  can  be  either  in  the  positive  or  ; 

\ 

j negative  direction,  the  sign  of  the  cross-range  i 

| ’ 

l error  (CRE)  must  be  preserved.  The  initial  missile 

| cross-range  error  can  be  in  either  direction;  ; 

I . consequently,  the  equation  becomes:  \ 

( • • 3 

i CRE  = CRE  jCRE|/2G  + CRE(O)  » 

E 

\ • * 

| This  equation  is  parabolic  as  depicted  in  Figure  \ 

I V-6 . 

J 

Also,  since  the  magnitude  of  the  missile  cross- 
range acceleration,  "G" , is  constant,  the  correspon- 
ding acceleration  vector,  "U" , is: 

U = ±G  + D = (-G) SIGN (error  function) 

i • • 

| , error  function  = EX  = CRE|CRE|/2G  + CRE 

This  drives  the  missile's  control  surfaces  in  the  proper 
direction,  a direction  which  is  opposite  to  the  sign  ;■ 

of  the  error  function,  thus  moving  the  missile  back 
toward  the  line  of  sight  with  a constant  magnitude  of 
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acceleration,  "G". 


A simplified  block  diagram  representing  the  command  to 
line-of-sight  guidance  system  is  presented  in  Figure  V-ll. 
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Figure  V-ll  Elementary  Command  to  Line-of-Sight  Block 
Diagram 


VI.  'BASIC'  COMMAND  TO  LINE-OP-SIGHT  SIMULATION 


A.  PROGRAM  DESCRIPTION 

Utilizing  the  concepts  previously  discussed,  this  program 
simulates  theoretically  optimum  surface-to-air  missile  maneuver 
against  a constant-velocity  target  in  two  dimensions.  The 
engagement  was  designed  to  occur  within  the  first  quadrant  of  a 
standard  Cartesian-coordinate  reference  system  with  the  ground 
tracker,  missile  lauching  unit,  located  at  the  origin.  The 
missile  was  assumed  to  be  "captured"  by  the  guidance  system 
approximately  one  second  after  missile  launch  or  at  a distance 
of  about  500  meters  from  the  ground  tracker  unit.  At  the  time 
of  missile  "capture"  an  initial  cross-range  error  of  ten  meters 
and  a cross-range  error  rate  of  50  meters  per  second  were  intro- 
duced to  demonstrate  the  ability  of  the  guidance  system  to  neu- 
tralize initial  errors.  It  should  be  noted  that  these  values 
of  error  were  selected  entirely  arbitrarily  and  that  they  could 
be  any  reasonable  values  assumed  to  exist  at  the  time  of  missile 
"capture" . 

The  target  was  "flown"  across  the  first  quadrant  from  a 
position  5000  meters  from  the  tracker  on  the  X-axis  (5000,0) 
toward  a position  5000  meters  from  the  tracker  on  the  Y-axis 
(0,5000) . 

In  order  to  use  a command  to  line-of-sight  guidance  scheme 
effectively,  the  missile  must  possess  a speed  advantage 
over  the  target.  Therefore,  in  this  model  the  missile  velocity 
was  chosen  to  be  500  meters  per  second  versus  a 250  meter-per- 
second  target.  At  sea  level  these  speeds  equate  to  approximately 


Mach  1.5  and  Mach  0.75  respectively.  The  missile,  also,  after 
being  launched  directly  toward  the  target,  was  maneuvered  using 
on-off,  thrust-vector  control  which  applied  a lateral  acceleration 
of  60  meters  per  second  squared  to  the  missile.  This  is  just 
greater  than  six  times  the  earth’s  gravitational  acceleration 
at  sea  level. 

In  general,  keeping  in  mind  the  scenario  above,  the  program 
operates  as  follows: 

The  equations  for  the  target's  movement  were 
written  and  appropriately  placed  into  the  program 
so  that  the  target  and  missile  angular  directions 
could  be  compared.  The  difference  between  these 
two  angles  creates  an  error  angle.  This  error 
angle  is  geometrically  converted  by  the  tracker 
computer  into  the  cross-range  error,  the  distance 
of  the  missile  from  its  intended  flight  path  - 
the  line  of  sight  to  the  target  in  this  case. 

That  error  distance  is  reduced  to  zero  by 
applying  a lateral,  thrust-vector-controlled 
acceleration  to  the  missile  in  the  correct 
direction  as  determined  and  commanded  by  the 

\ ground  tracker  unit.  The  optimum  parabolic 

I 

l 

^ switching  boundary,  discussed  in  Section  V,  was 

used  in  the  tracker  computations.  In  order  to 
determine  this  cross-range  error,  the  tracker 
estimates  the  missile's  range  by  the  elapsed  time 
of  flight  and  the  missile's  velocity  profile. 

Also,  to  simplify  the  program  somewhat,  the 
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missile's  velocity  was  assumed  to  be  constant 


in  a radial  direction  from  the  tracker.  When  the 


cross-range  error  is  reduced  to  an  acceptable  dis- 


tance, the  amount  of  thrust-vector  control  is 


limited  to  prevent  rapid  changes  between  full 


thrust-vector  deflections  in  both  directions.  If 


the  error,  again,  becomes  unacceptably  large,  the 


full  deflection  positions  are  utilized.  Missile 


position  is  determined  by  state  variable  analysis 


using  the  INTEG2  subroutine  for  integration  of  the 


lateral  thrust-vector-controlled  missile  acceleration. 


The  use  of  this  subroutine  is  explained  in  Appendix 


A.  After  two  integrations,  the  missile's  new  cross' 


range  error  is  determined  and,  then,  this  distance  is 


used  with  estimated  missile  range  to  provide  a revised 


missile  angular  position.  This  process  is  repeated 


in  increments  of  five  thousandth  of  a second  from 


the  time  of  missile  "capture"  by  the  tracker,  a 


time  of  approximately  one  second,  until  maximum 


missile  range  is  attained  after  ten  seconds  of 


flight.  Thus,  the  missile  is  driven  continually 


to  a position  upon  the  line  of  sight  to  the  target 


The  program  variables  are  "frozen"  when  the  missile 


reaches  the  closest  point  of  approach  to  the  target 


if  the  value  of  constant  C(l)  is  one.  If  this 


value  is  zero,  the  program  computes  until  the  time 


of  missile  flight  reaches  ten  seconds. 
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The  program  was  written  in  the  FORTRAN  IV  language.  A 
detailed,  annotated  copy  of  the  program,  including  a table  of 
program  variable  definitions,  is  provided  following  the 
narrative  portion  of  this  paper. 
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Figure  VI-1  'Basic'  Computer  Simulation  Block  Diagram 

B.  SYSTEM  BLOCK  DIAGRAM 

The  command  to  line-of-sight  computer  simulation  block 
diagram.  Figure  VI-1,  shows  how  the  elementary  block  diagram 
of  this  guidance  technique,  Figure  V-ll,  was  expanded  upon  in 
order  to  develop  a computer  model  for  the  system.  As  will  be 
seen,  while  appearing  to  be  quite  simple  in  concept,  the 
simulation  becomes  increasingly  more  complicated  as  its  sophis- 
tication is  enhanced. 

C.  SIMULATION  RESULTS 

The  plotted  output  of  the  phase  plane,  Figure  VI-2,  shows 
how  the  cross-range  error  and  cross-range  error  rate  are  driven 
to  a near  zero  condition  along  the  optimum,  parabolic,  switching 
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boundary.  Once  the  cross-range  error  is  less  than  one  meter, 


the  magnitude  of  the  thrust  vector  control  is  reduced  by  a 


sixth  to  minimize  the  possibility  of  over-controlling  the 


missile. 
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Figure  VI-2  'Basic'  Simulation  Phase  Plane 


Figure  VI-3  presents  the  "X"  and  "Y"  positions  for  the 


missile  and  the  target  versus  time.  As  can  be  seen  the  inter- 


cept occurs  in  just  under  eight  seconds  after  missile  launch. 


Figure  VI-4  indicates  the  relative  positions  of  the 


missile  and  target  with  the  terminal  points  for  each  flight 


profile  stopped  at  the  closest  point  of  approach.  An  intercept 


did  occur. 


In  order  to  ensure  that  the  missile's  cross-range  error 


from  the  line  of  sight  was  damped,  Figure  VI-5  was  produced  to 


show  this  cross-range  error  versus  time  of  flight.  The  initial 


errors,  applied  at  the  time  of  missile  "capture”  by  the 
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guidance  system  (corresponding  to  a time  of  approximately  one 
second  in  the  figure),  were  reduced  sufficiently  in  less  than 
two  seconds. 


D.  EFFECTS  OF  DECEPTIVE  JAMMING 

A readily  apparent  deceptive  jamming  technique,  which  might 
be  effective  against  this  class  of  guidance  system,  would  be  to 
employ,  aboard  the  target  aircraft,  an  emitter  which  imitates 
the  signal  that  is  used  by  the  tracker  to  locate  the  missile’s 
angle  from  the  reference  plane.  This  could  be,  for  example, 
depending  upon  how  the  missile  is  tracked  in  angle,  an  imitation 
of  the  missile's  transponder  radio- frequency  signal  or  a dup- 
lication of  the  missile's  infrared  emissions.  In  any  case  the 
deceptive  jamming  signal  from  the  target  aircraft  should  be  of 
sufficient  strength  to  ensure  that  the  missile  tracking  system 
within  the  ground- located  tracker  considers  the  target  to  be 
the  missile.  There  are,  most  assuredly,  counter-countermeasures 
which  could  be  applied  to  negate  this  deceptive  technique;  how- 
ever, this,  currently,  would  be  one  method  of  generating  a zero 
error  signal  between  the  missile  and  target  angles.  Additionally, 
this  deceptive  countermeasure  provides  a method  to  demonstrate 
how  the  computer  simulation  drives  the  missile  with  inputs  of 
erroneous  guidance  information. 

In  order  to  implement  this  deceptive  jamming  technique,  it 
was  assumed  that  the  target  aircraft  jamming  equipment  became 
effective  when  the  missile  was  approximately  half  the  distance 
to  the  target.  This  would  equate  to  about  four  seconds  of 
missile  flight  time.  Therefore,  since  this  jamming 
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countermeasure  system  simply  zeros  the  error  angle,  orly  one 
additional  equation  was  necessary  to  modify  the  previously 
discussed  "basic"  command  to  line-of-sight  guidance  model.  The 
equation,  IF (T.GE.4.0/E=0. , was  placed  into  the  program  following 
the  computational  step  for  the  calculation  of  the  angular  error, 
E=SIGM-SIGT. 

The  resulting  missile  trajectory  and  other  parameters  are 
readily  apparent  in  Figure  VI-6  through  VI-8.  As  shown  in 
Figure  VI-6,  the  error  was  reduced  to  zero  by  the  time  that 
jamming  commenced.  After  deceptive  jamming  was  initiated,  the 
cross-range  error  rapidly  increased  causing  a miss  of  over  300 
meters.  Needless  to  say,  because  of  the  short  times  of  flight 
involved,  the  countermeasures  must  be  employed  in  a timely 
manner.  Otherwise,  the  missile  could  approach,  easily,  within 
lethal  range  of  the  target  aircraft. 
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Figure  VI-6  ’Basic*  Simulation  Phase  Plane  When 
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Figure  VI- 7 ’Basic'  Simulation  Missile  & Target 
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Figure  VI-8  'Basic'  Simulation  Missile  Cross-Range  Error 
When  System  Deceptively  Jammed 
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VII.  ’LEAD-ANGLE'  COMMAND  TO  LINE-OF-SIGHT  SIMULATION 


A.  PROGRAM  DESCRIPTION 


This  'lead-angle'  simulation  uses  the  same  target  profile  as 


that  of  the  'basic'  command  to  line-of-sight  simulation.  This 


allows  a comparison  of  the  two  simulations  to  a certain  degree. 


One  parameter,  which  can  not  be  compared  accurately,  is  the  time 


to  impact.  This  occurs  because  the  radial  missile  velocity  from 


the  tracker  was  assumed  to  be  a constant  value  - in  this  case 


500  meters  per  second.  Therefore,  intercepts  will  take  place 


at  about  the  same  time  regardless  of  missile  maneuvers.  A 


general  comparison  can  be  made,  however,  by  mentally  comparing 


the  plotted  missile  flighr  profiles,  the  longer  profile  taking 


the  most  time  in  an  actual  engagement.  Also,  in  order  to  avoid 


redundancy,  discussion  will  be  limited  to  modifications  of  the 


previously  described  computer  model.  Only  a few  changes  were 


required  to  produce  a 'lead-angle'  variant  for  the  command  to 


line-of-sight  guidance  system. 


In  this  guidance  scheme,  the  spacial  positions  of  both  the 


target  and  the  missile  must  be  accurately  and  continuously 


known  in  order  to  compute  an  advanced  impact  position.  This  is 


accomplished  by  calculating  the  distance  from  the  missile  to 


the  target  and  their  closing  velocity.  From  this  data,  a time 


remaining  until  impact  can  be  estimated.  This  remaining  time 


of  flight,  then,  is  combined  with  the  target's  current  velocity 


vector  and  position  to  ascertain  the  expected  impact  point.  A 


synthetic  line  of  sight  toward  this  location  is  generated  by 


the  ground  tracker  unit  computer  and  the  missile  is  commanded 
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to  "fly"  up  this  flight  path  instead  of  along  the  line  of  sight 
directly  toward  the  target.  The  guidance  scheme  used  to  main- 
tain the  missile  on  course,  the  synthetic  line  of  sight,  is, 
otherwise,  identical  to  that  of  the  "basic"  simulation  previously 
described.  Once  the  missile  passes  the  closest  point  of 
approach  to  the  target,  this  simulation  ceases  to  be  valid.  This 
is  quite  all  right  because,  once  the  missile  has  passed  the 
target,  no  maneuver  would  be  sufficient  to  effectively  engage 
the  aircraft. 


B.  SYSTEM  BLOCK  DIAGRAM 

The  "lead-angle"  command  to  line  of  sight  computer  simulation 
block  diagram  is  contained  in  Figure  VII-1.  It  is  essentially 
a modification  of  the  "basic"  simulation  block  diagram,  Figure 
VI-1.  As  described  previously,  target  and  missile  positions 
plus  their  velocities  must  be  known  for  the  guidance  scheme  to 
be  viable. 
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Figure  ^II-l  'Lead-Angle*  Computer  Simulation  Block 
Diagram 
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C.  SIMULATION  RESULTS 


This  simulation  supports  the  idea  that  the  use  of  a "lead- 
angle"  variant  to  the  "basic"  command  to  line-of-sight  guidance 
system  enhances  the  effectiveness  of  the  elementary  guidance 
technique*  As  can  be  seen  in  Figures  VII-2  through  VII-5,  the 
missile  was  "captured"  by  the  guidance  system  at  a time  of 
approximately  one  second  from  launch.  The  initial  cross-range 
and  cross-range-rate  errors  were  quickly  reduced  and  the  missile 
"flew"  out  the  synthetic  line-of-sight  to  successfully  intercept 
the  target  aircraft.  Of  course,  if  the  target  employs  electronic 
countermeasures  such  as  denial  jamming  of  the  target  tracking 
equipment,  this  "lead-angle"  method  would  prove  to  be  useless 
because  accurate  positioning  information  must  be  maintained  to 
compute  the  estimated  impact  point.  Also,  if  the  target 
maneuvers  sufficiently  to  cause  the  impact  point  to  vary  at 
high  rates,  the  synthetic  line  of  sight  may  shift  at  too  high 
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a rate  for  the  missile  to  follow.  Thus,  the  missile  may  not 
be  able  to  physically  engage  the  impact  point  but  might  be 


able  to  still  effectively  attack  the  target  by  the  "basic" 
guidance  scheme.  Consequently,  in  order  to  capitalize  upon 
the  advantage  of  the  "lead-angle"  variant  (an  increased  effee 
tive  missile  envelope  against  a relatively  slow  maneuvering 
target  which  is  not  employing  electronic  countermeasures) , 
it  seemed  that  a combined  guidance  technique  should  be  con- 
sidered. Therefore,  the  "consolidated"  simulation  discussed 
in  the  next  section  was  contrived. 
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Figure  VII-5  'Lead-Angle'  Simulation  Missile  Cross-Range 
Error 
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VIII.  ’CONSOLIDATED*  COMMAND  TO  LINE-CF-SXGHT  SIMULATION 


A.  PROGRAM  DESCRIPTION 

When  the  two  previous  simulations  were  combined,  a major 
problem  surfaced.  That  problem  was  basing  the  missile's  new 
position  upon  two  different  references,  the  line  of  sight  to 
the  target  in  the  case  of  the  "basic”  model  and  the  synthetic 
line  of  sight  to  the  impact  point  in  the  case  of  the  "lead-angle" 
variant.  This  incompatibility  caused  the  missile's  position 
to  jump  radically  as  the  "consolidated"  simulation  attempted 
to  shift  alternately  between  the  two  reference  lines  as  the 
jamming  situation  varied.  Therefore,  a single  reference  had 
to  be  implemented  for  continuous  use  regardless  of  whether 
jamming  was  or  was  not  encountered.  The  X-axis  of  the  Cartesian- 
coordinate  system  was  selected  since  various  angular  measurements 
were  based  upon  it  already.  Also,  this  has  proved  to  provide 
a more  accurate  model  of  the  missile’s  optimum  flight  profile 
because  the  reference  does  not  move  during  the  time  of  missile 
flight.  If  a strictly  "basic"  simulation  should  be  desired, 
jamming  is  considered  to  be  effective  during  the  entire  engage- 
ment. Conversely,  the  absence  of  jamming  during  the  missile's 
flight  produces  a total  "lead-angle"  intercept  of  the  target 
aircraft.  Of  course,  different  duty  cycles  for  jamming  can  be 
implemented  by  the  target  during  a single  computer  run  to  test 
the  ability  of  the  missile  system  to  transition  between  the  two 
attack  modes.  A limited  number  of  different  jamming  combinations 
were  tested  and  are  discussed  later. 
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In  order  to  implement  jamming,  the  jamming  time  increments 
are  read  into  the  simulation  through  five  time  intervals 
specified  on  the  first  five  data  cards.  These  cards  preceed 
those  required  by  the  INTEG2  subroutine  as  explained  in 
Appendix  A.  The  time  each  jamming  interval  begins  is  punched 
into  columns  1-10  of  each  card  in  F-format.  The  times  that 
jamming  ceases  are  entered  in  columns  11-20.  Should  a no- 
jamming condition  be  desired,  then,  all  five  cards  are  left 
blank  but  they  must  not  be  ommitted.  Total  jamming  can  be 
obtained  by  making  one  card  cover  the  entire  time  of  flight. 
For  example,  the  first  card  could  have  columns  1-10  blank  and 
the  value  11.0  entered  into  columns  11-20.  This  produces 
jamming  from  time  zero  until  a time  of  eleven  seconds,  which 
is  well  after  the  end  of  ■'die  engagement.  The  remaining  four 
cards  are  blank  and  must  be  inserted. 

This  simulation  "drives"  the  missile  in  the  "Y"  direction 
with  sufficient  "Y"  velocity  to  reduce  the  cross-range  error 
to  zero.  The  missile's  X-direction  velocity  component  is  used 


to  maintain  a constant,  outward,  radial  velocity,  "V  ",  toward 

m 


the  target.  In  this  simulation,  as  in  its  predecessors,  that 


radial  velocity,  "V  ",  is  assumed  to  be  500  meters  per  second. 


Therefore,  the  velocity  of  the  missile  in  the  "X"  direction  is 
related  to  that  obtained  from  the  forcing  function  in  the  "Y" 
direction.  This  is  shown  in  Figure  VIII-1. 


Figure  VIII-1  'Consolidated'  Simulation  Missile  Velocity 
Vectors 
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The  missile's  Y-direction  velocity  is  obtained  from  the  state 
variable  analysis  through  a single  integration  of  the  forcing 
function,  "U" , which  is  the  missile's  Y-direction  acceleration. 
This  acceleration  is  the  "Y"  component  of  the  cross-range  error 
acceleration  achieved  by  thrust-vector  control  of  the  missile  as 
shown  in  Figure  VIII-2. 
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Figure  VIIX-2  'Consolidated'  Simulation  Missile 
Acceleration  Vectors 


The  missile's  Y-coordinate  position,  then,  is  the  second  inte- 
gration of  the  forcing  function.  From  these  results,  the  X- 
direction  components  were  determined  since  the  velocities  must 
yield  a constant,  radial  velocity  of  500  meters  per  second.  In 
order  to  ensure  a correct  integration  process,  the  missile's 
"Y"  velocity  was  not  limited  to  500  meters  per  second;  however, 
it  never  exceeded  that  value  in  any  of  the  computer  runs  tested. 
It  should  be  noted  that  the  program  will  zero  the  acceleration 
and  the  velocity  of  the  missile  in  the  "X"  direction  if  excessive 
missile  "Y"  velocity  is  experienced.  Also,  the  results,  then, 
will  become  invalid;  consequently,  the  velocity  of  the  missile. 
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«V , should  be  increased  or  the  target 1 s speed  reduced  to 
obtain  a theoretically-possible  engagement  scenario. 

The  targets,  in  the  previous  simulations,  "flew"  a constant- 
velocity,  linear  course.  In  this  simulation  the  target  was 
maneuvered  sinusoidally  in  the  "Y"  direction  for  the  purpose 
of  observing  the  effects  of  a moving  impact  point.  The  maximum 
velocity  of  the  target  in  the  "Y"  direction  is  64  meters  per 
second  with  a half  period  of  four  seconds.  The  target's  "X" 
velocity  was  a constant,  negative  210  meters  per  second.  This 
did  produce  some  interesting  results.  First,  this  type  of  tar- 
get maneuver  caused  the  impact  position  to  move  at  even  higher 
rates  than  the  target  moved.  Thus,  in  some  cases,  the  missile 
may  not  be  able  to  engage,  adequately,  the  impact  point  which 
is  moving  at  a very  high  velocity.  It,  however,  may  be  able  to 
attack  the  target  itself  by  utilizing  the  "basic"  command  to 
line-of-sight  approach.  To  determine  which  situation  exists 
necessitated  the  use  of  the  impact  point's  velocity  to  compute 
the  synthetic  line  of  sight,  the  angle  upon  which  the  missile 
"flies"  in  a non- jamming  environment.  If  this  impact-point 
velocity  exceeded  eight- tenths  of  the  missile's  velocity,  then, 
the  missile  was  commanded  to  "fly"  directly  toward  the  target's 
present  position  by  the  "basic"  guidance  method.  This  was  in 
lieu  of  utilizing  the  "lead-angle"  approach  with  only  an  un- 
stable impact  position  and,  thus,  an  erratic  or  too  quickly 
moving  synthetic  line  of  sight. 

In  order  to  provide  a longer  time  of  flight,  the  target's 
initial  position  was  moved  to  coordinates  (6000,1000)  and  the 
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target  was  "flown"  inbound  at  210  meters  per  second  while 
sinusoidally  maneuvering  about  the  1000  meter  Y-ordinate.  Also, 
in  regard  to  the  short  time  of  flight  involved,  excessively 
large  initial  errors  in  either  distance  or  rate  should  be 
avoided  because  of  the  time  required  to  neutralize  these  errors. 
For  example,  since  the  missile's  maximum  acceleration  is  only 
150  meters  per  second  squared,  initial  erroneous  conditions  of 
250  meters  and  250  meters  per  second  for  the  missile's  Y-direction 
components  require  many  seconds  to  correct.  However,  with 
initial  conditions  of  150  meters  and  150  meters  per  second,  the 
simulation  proved  the  missile  to  be  effective  against  this 
target  over  a full  scope  of  jamming  duty  cycle  variations  - 
from  none  to  continuous  denial  jamming. 

The  initial  angle  to  the  missile  at  the  time  of  "capture" 
by  the  guidance  system,  a time  of  one  second  or  500  meters  from 
the  launcher-tracking  unit,  is  determined  by  the  initial  con- 
ditions, X(l)  and  X (2) . As  can  be  recalled,  X(l)  is  the 
missile's  "Y"  position  at  the  time  that  the  weapon  is  500  meters 
radially  from  the  launcher.  Therefore,  the  angle  is  the  inverse 
sine  of  X(l)  divided  by  500  meters.  In  order  to  attain  this  "Y" 
position  in  one  second,  the  missile's  "Y"  velocity,  then,  was 
approximately  the  "Y"  distance  per  second.  Consequently,  in 
this  simulation,  initial  conditions  of  100  meters  for  X(l)  and 
100  meters  per  second  for  X(2)  equate  to  an  angle  of  11.54 
degrees  above  reference. 

Also,  the  magnitude  of  the  missile's  maximum  lateral 


acceleration,  "G",  was  increased  to  150  meters  per  second 
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squared.  This  provided  a sufficient  forcing  function  to  allow 
engagement  of  the  target.  This  was  possible  even  with  increased 
initial  errors  at  the  time  of  missile  "capture"  by  the  guidance 
system;  with  a maneuvering  target  accelerating  at  just  over  50 
meters  per  second  squared,  at  maximum,  in  the  "Y"  direction; 
and  with  various  jamming  duty  cycles  from  zero  to  one.  As  in 
the  previous  simulations,  the  thrust- vector  control  of  the 
missile  was  reduced  by  one-sixth;  however,  this  was  done  only 
when  the  sum  of  the  absolute  values  of  cross-range  error  and 
error  rate  was  less  than  or  equal  to  four  vice  the  value  of  one, 
which  was  used  previously.  The  error  tolerance  was  increased  to 
provide  a larger  "idling"  area  for  the  missile's  thrust-vector 
control  because  of  the  increased  forcing  function  utilized. 

The  actual  computer  program  for  the  simulation,  with 
comments  included,  follows  the  appendices. 
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B.  SYSTEM  BLOCK  DIAGRAM 

The  block  diagram  of  the  "consolidated"  simulation  is  con- 
tained in  Figure  VIII-3 . Although  it  appears  formidable  at 
first  glance,  the  diagram  merely  combines  the  concepts  of  the 
two  previous  simulations  with  slight  modifications  to  account 
for  the  use  of  a new  reference  from  which  to  compute  the 
missile's  location.  If  the  tracking  unit  ascertains  that  it  is 
contending  with  either  denial  jamming  or  a highly-maneuvering 
target,  then,  only  the  upper  portion  of  the  tracker  diagram  is 
used  for  guidance.  Conversely,  if  both  a non- jamming  condition 
and  a relatively  slowly-maneuvering  target  situation  are 
experienced,  the  lower  portion  of  the  tracker  diagram  produces 
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guidance  for  the  missile. 

The  simulation  does  produce  necessary  calculations  in  all 
areas  of  the  tracker  diagram  continually  throughout  the  engage- 
ment, but  only  generates  commands  to  the  missile  from  the 
appropriate  section  of  the  model. 

After  the  program  reads  the  data  cards,  the  times  of 
jamming  are  known  and  the  appropriate  control  function  within 
the  tracker's  computer  is  implemented  to  control  the  missile. 

The  magnitude  of  the  impact  point's  velocity  is  considered,  also, 
to  ascertain  the  stability  of  the  synthetic  line  of  sight  and, 
thus,  to  determine  whether  target  maneuvers  make  the  ''lead- 
angle''  approach  of  guidance  inappropriate. 


C.  SIMULATION  RESULTS 

| Three  engagements  were  computed  under  different  jamming 

conditions.  They  were  continuous  jamming,  intermittent  jamming  II 

having  an  approximate  duty  cycle  of  one-half,  and  no  jamming. 

During  intermittent  jamming,  the  ground-located  tracker  unit 
was  denied  range  information  within  the  following  time  increments 
after  missile  launch:  1. 0-3.0,  4. 0-5.0,  6. 6-7.0,  7. 4-7. 8,  and 
8. 2-8. 6 seconds. 

Utilizing  the  scenario  which  was  previously  discussed,  all 
three  attacks  resulted  in  successful  intercepts  of  the  target. 

As  can  be  observed  in  Figure  VIII-4,  the  best  flight  profile 
of  the  three  simulated  runs  occurred  when  total  jamming  was 
experienced.  This  jamming  caused  the  missile  to  "fly" 
directly  toward  the  target. 
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Figure  VIII-4  'Consolidated'  Simulation  Mssile  & Target 
Flight  Paths  When  Continue1  O.y  Jammed 
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When  no  jamming  was  encountered,  the  flight  profile  became 


more  exaggerated  because  the  target  was  maneuvering,  thus 


causing  the  estimated  impact  point  to  move  at  an  even  higher 


rate.  Recall  that,  in  a non- jamming  environment,  the  missile 


will  "fly"  for  the  impact  point  vice  the  target.  Figure  VIII-5 


shows  this  profile.  Finally,  with  intermittent  jamming,  the 


missile's  flight  profile,  Figure  VIII-6,  became  more  erratic 


but  did  not  reach  the  more  extreme  deviations  experienced  in  the 


un jammed  case.  Therefore,  it  appears  that  the  best  method  would 


be  to  use  the  "basic"  command  to  line-of-sight  guidance  technique 


against  a maneuvering  target  regardless  of  the  jamming  condition 


in  effect.  Also,  as  can  be  observed  in  Figure  VIII-6,  the  flight 


profile  during  intermittent  jamming,  the  amount  of  erratic 


variation  in  the  missile's  movement  decreased  considerably  as  it 


closed  the  target.  This  occurred  because  the  difference  between 


the  synthetic  line  of  sight  and  the  line  of  sight  becomes 


smaller  with  time  until  they  are  equal  in  value  at  the  time  of 


impact.  Consequently,  the  "lead-angle"  approach  should  be 


considered  for  utilization  against  a non-maneuvering  target 


when  adequate  range  information  is  available,  i.e.,  when  denial 


jamming  is  not  encountered. 


In  order  to  reduce  the  amount  of  movement  in  the  missile's 


thrust-vector  controllers  and,  thus,  to  ensure  that  they  do 


not  "beat  themselves  to  death"  mechanically  while  trying  to 


reduce  a trivial  cross-range  error,  a reduced  amount  of 


deflection  was  applied  when  the  cross-range  error  was  relatively 


small.  This,  then,  was  a trade-off  of  allowing  some  increase 
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in  cross-range  error  for  improved  thrust-vector  control 
mechanical  reliability.  Figure  VIII-7  shows  how  the  crosr- 
range  error  was  reduced  without  using  a smaller  amount  of  thrust 
vector  control  deflection.  Keep  in  mind  that  the  control  sur- 
faces in  the  missile  were  shifting  rapidly  from  a total  positive 
to  a total  negative  direction.  In  Figure  VIII-8,  a cross- 
range error  of  four  meters  was  accepted  as  a reasonable  maximum 
miss  distance  and  the  amount  of  control  was  decreased  when  the 
error  was  within  chat  value.  This  assured  improved  mechanical 
reliability  of  the  system.  As  can  be  seen  in  Figure  VIII-8, 
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Figure  VIII-7  Figure  VIII-8 

'Consolidated*  Simulation  'Consolidated*  Simulation 

Missile  Cross-Range  Error  Missile  Cross-Range  Error 

With  Full  Missile  Control  With  Reductions  in  Missile 
Surface  Deflections  Control  Surface  Deflections 


at  the  time  of  impact  the  cross-range  error  was  in  the  process 
of  being  reduced  from  a value  of  about  four  meters.  The  phase 


planes.  Figures  VIII-9  and  VIII-10,  also,  show  the  effects  of 


this  trade-off.  Figure  VIII-9  was  the  result  obtained  when 


using  no  reduction  in  the  amount  of  thrust-vector  control; 


while  Figure  VIII-10  shows  the  effect  of  allowing  a maximum 


four-meter  error  with  a reduced  amount  of  lateral  thrust. 


IX.  CONCLUSIONS 


These  simulations  support  the  premise  that  consideration 
should  be  given  to  utilizing  a "lead-angle"  variant  in  con- 
junction with  the  "basic"  command  to  line-of -sight  guidance 
technique.  This  would  increase  the  missile  system's  effective 
engagement  envelope  against  a non-maneuvering  target  that 
employs  low-duty-cycle  denial  jamming.  If  the  target  utilizes 
higher  duty  cycle  jamming  or  begins  to  make  evasive  maneuvers, 
then,  the  system  can  revert,  when  necessary,  to  the  "basic" 
guidance  mode  even  after  the  missile  has  been  launched. 

This  new  approach  appears  to  be  generically  applicable 
to  missiles  of  the  command  to  line-of-sight  guidance  class. 

If  a range- tracking  capability  already  is  contained  within 
the  system,  it  should  be  possible  to  implement  this  improved 
guidance  technique  through  modifications  to  the  existing 
guidance  computer. 
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APPENDIX  A 


DESCRIPTION  AND  USE  OF  SUBROUTINE  INTEG2 


INTEG2  provides  a simple  means  for  a Runga-Kutta  solution 
of  ordinary  differential  equations.  Input  and  output  routines 
(including  graphical  point  plots  on  the  versatec  plotter)  are 
built  in. 

The  subroutine  has  been  compiled  and  added  to  the  FORTRAN 
library,  SYS1.MPSLIB.  The  program  is  flexible  and  self  con- 
tained. It  will  seldom  be  desirable  to  build  the  subroutine 
into  a more  complex  program. 

It  is  assumed  that  the  equations  to  be  solved  have  been 
reduced  to  'N'  simultaneous  first  order  differential  equations 
of  the  form: 

DX(I)/DT  » F(X(1) ,X(2) ,. . . ,X(N) ,T) , I = 1,2, . . . ,N 
where  the  function  'Fr  may  be  nonlinear.  The  user  must  trans- 
late these  'N'  equations  into  FORTRAN  as  in  the  following  second- 
order  example: 

XD0T(1)  = X (2) 

XDOT ( 2 ) = C (1) *X(1)-0.2*X(2)+C(2) *SIN(T*T) 

C(l)  and  C(2)  are  constants  which  are  entered  on  data  cards. 

To  solve  these  equations  the  OS/360  deck  must  be  of  the 
following  form: 


// STANDARD  JOB  CARD  (GREEN) 

//  EXEC  FORTCLGV, REGION. G0=150K 
//FORT. SYS IN  DD  * 

DIMENSION  X(30) ,XDOT(30) ,C(15) 


% . >1 

r ' 


I?'  A 


if:-  i 


! ! 

*.  \ 


1 


C(10)  = 1. 


1 CALL  INTEG2 (T,X,XDOT,C) 


USERS  FORTRAN  EQUATIONS 


GO  TO  1 


//GO. SYS IN  DD  * 


DATA  CARDS  AS  DESCRIBED  BELOW 


/*  (ORANGE) 


The  number  of  XDOT  equations  is  fN',  which  must  not  exceed 


30.  These  equations  are  supplied  by  the  user,  who  defines 


each  'XDOT'  in  terms  of  the  dependent  variables  X(l)  through 


X (N)  and  the  independent  variable  'T'.  In  writing  these 


equations  the  user  may  introduce  at  his  convenience: 


any  unsubscripted  variables; 


the  subscripted  variables  X(I) , where  the  variable 


I'  must  be  less  than  or  equal  to  30; 


the  constants  C(I)  to  be  entered  as  data,  where  the 


variable  ' I*  must  be  less  than  or  equal  to  eight; 


any  normal  FORTRAN  technique  or  function;  and 


routines  from  any  source  library  or  user-supplied 


subroutines . 


Note  1:  The  use  of  'auxiliary'  X(I),  when  the 


variable  ' I ' is  greater  than  the  number  of  equations 


to  be  integrated,  does  not  alter  the  value  of  'N' , 


62 


the  order  of  the  equations. 

Note  2:  Loops,  either  with  or  without  a 'DO'  state- 
ment, are  best  avoided. 

Note  3:  'IF'  statements,  provided  that  they  do  not 
create  a loop,  can  be  used  to  transfer  control 
within  the  user’s  equations.  For  example,  the 
statement,  IF(T.GT.10.)C(3)=0. , would  cause  C(3) 
to  take  the  value  zero  for  all  'T'  greater  than  ten. 

The  constants  C(l)  through  C(8)  may  be  used  as  described  in 
the  above  examples,  and  are  read  in  from  a data  card.  C(10) 
must  never  be  used,  except  as  indicated  in  the  standard  OS/360 
deck  above.  C(ll)  and  C(12)  control  the  output.  For  example, 
if  the  statements,  C(ll)  = 10.0  and  C(12)  * 2.0,  are  added  to 
the  users  FORTRAN  equations,  every  tenth  integration  step  will 
be  printed  out,  and  every  second  step  will  be  plotted.  (If 
not  set  by  the  user,  default  values  of  20  and  five  apply.) 

C (13)  can  be  similarly  used  to  modify  the  step  size  of  the 
numerical  integration  which  is  more  usually  defined  on  a data 
card. 

Only  the  independent  variable  'T*  and  the  variables  X(l) 
through  X(30)  can  be  output.  Therefore,  to  output  a quantity 
which  is  not  one  of  the  ’N’  dependent  variables.  The  user 
must  introduce  an  'auxiliary'  variable  X(I),  where  the  variable 
'I'  must  be  less  than  31  but  greater  than  'N'  - for  example, 
by  introducing  the  FORTRAN  equation,  X(27)  = X(3)*X(3),  the 
square  of  X(3)  can  be  output  by  the  program  as  X(27) . Normally, 
a line  of  printout  is  generated  after  every  20  integration 
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steps,  and  one  graph  point  (per  curve)  after  every  five  inte- 
gration steps.  These  values  can  be  modified  as  indicated  in  the 
previous  section.  A run  is  terminated  if  the  run  exceeds  450 
lines  or  if  more  than  900  graph  points  are  generated  per  curve. 

The  accuracy  of  the  numerical  integration  is  controlled  by 
the  choice  of  step  size.  Unfortunately  there  is  no  simple 
method  for  making  this  choice  other  than  by  trial  and  error. 

As  a first  guess,  a step  size  of  1/1000  of  the  total  range  of 
interest  of  the  independent  variable  'T'  is  usually  reasonable. 
If  one  or  more  of  the  variables  change  rapidly  in  part  of  the 
range  of  ' T',  then  a finer  step  size  may  be  necessary  for  all 
or  part  of  the  run.  The  only  way  to  be  reasonably  sure  of  the 
solution  accuracy  is  to  re-run  the  problem  with  a smaller  step 
size  and,  then,  to  compare  the  results.  A maximum  of  4500 
integration  steps  is  permitted  in  any  one  run. 

Data  cards  are  assembled  as  follows: 


First: 


Second : 


Third: 


The  user's  job  identification  label  is 
punched  in  columns  1-42. 

The  number  of  runs  to  be  processed,  which 
must  be  less  than  ten,  is  punched  in  column 
one.  The  run  number  together  with  the  user's 
job  identification  label  is  placed  on  all 
output. 

The  order  of  the  differential  equations, 
which  must  be  less  than  31,  is  punched 
with  two  digits  in  columns  one  and  two  - 
for  example,  '03'. 


Fourth:  Initial  and  final  values  of  the  independent 

variable  (TI  and  TF)  and  the  step  size  (DT) 
appear  in  the  order  TI  - DT  - TF  and  are 
punched,  with  decimal  points,  in  columns 
1-10,  11-20  and  21-30.  Also,  it  is  possible 
to  process  the  integration  in  either  two  or 
three  segments  of  differing  step  size.  The 
data,  in  columns  of  ten  as  above,  then  takes 
the  form  TI  - DTI  - TF1  - DT2  - TF  or  TI  - 
DTI  - TF1  - DT2  - TF2  -DT3  - TF . 

Fifth:  The  values  of  the  constants  C(l),  C(2) 

C(8)  are  punched,  with  decimal  points,  in 
columns  1-10,  11-20 ,... ,71-80 . Ten  blank 
columns  may  be  used  for  a zero  or  unused 
constant.  If  no  constants  are  used,  this 

entire  card  will  be  blank  but  it  must 

not  be  omitted. 

Sixth:  The  initial  values  of  X(l) , X ( 2) , . . . ,X (N) 

are  punched,  with  decimal  points,  in 
columns  1-10,11-20,...  . Additional  cards 
are  required  if  'N'  is  greater  than  eight. 

Ten  blank  columns  may  be  used  for  a zero 
initial  value. 

Seventh:  This  card  controls  the  choice  of  variables 

for  printout.  Each  group  of  ten  columns, 
1-10,11-20, .. .up  to  71-80,  may  be  used  to 
specify  a column  heading  (eight  characters) 


'-•v-Vc  r,xi ?•**■**<&  a* 


65 


SSSSti 


J a«MM*ii«»^''  " * "*  ? ' * '"fS  ^ - -C 


Eighth : 


and  a two-digit/  right- justified  subscript 
identifying  the  corresponding  variable  - 
for  example:  TIME  00  and  SPEED  03  in 
columns  1-10  and  11-20  would  cause  the 
independent  variable  1 T ' / represented  on 
this  card  by  the  subscript  00,  to  be 
printed  out  under  the  column  heading  'TIME' 
and  the  variable  X(3)  to  be  printed  out 
under  the  heading  'SPEED'.  No  more  than 
eight  variables  can  be  output  during  one 
run.  If  no  printout  is  desired,  this  card 
must  be  blank. 

This  card  controls  the  choice  of  variables 
for  graph  output.  Up  to  four  curves  can  be 
plotted,  either  all  on  separate  graphs  or 
all  on  a single  graph.  Each  group  of  20 
columns  1-20,  21-40 ,.. .specifies  a curve  - 
for  example:  SPEED  VS.  TIME  0300  (in  columns 
1-20)  would  cause  X(3)  to  be  plotted 
vertically  against  the  independent  variable 
'T'  horizontally.  (Again,  the  subscript  00 
represents  the  independent  variable  'T'.) 

The  graph  output  would  be  labeled  'SPEED  VS. 
TIME'.  Note  that  the  first  16  columns  of  each 
group  are  used  for  the  label,  the  seventeenth 
and  eighteenth  for  the  'Y'  ordinates,  and  the 
last  two  for  the  'X'  ordinates. 


w 


V - 


If  the  curves  are  to  be  on  separate  graphs,  each  graph 
must  have  a label.  If  the  curves  are  to  be  all  on  one  graph, 
only  the  first  label,  columns  1-16,  must  be  provided.  The 
other  labels  must  be  blank.  The  entire  card  must  be  blank  if 
no  graphical  output  is  required. 

If  several  solutions  are  required  of  the  same  equations 
but  with  different  constants  or  initial  conditions,  the  number 
of  runs,  specified  on  the  second  data  card,  can  be  increased. 

For  all  runs  after  the  first,  only  the  fourth  through  the  last 
data  cards  are  to  be  supplied.  Except  for  the  data  on  the 
first  three  cards,  and  the  equations  themselves,  no  information 
is  retained  between  runs. 

Caution  should  be  exercised  by  the  user  because  too  large 
an  integration  step  size  can  result  in  unacceptable  errors  in 
the  solution,  and  even  in  instability  in  the  integration  process. 

This  subroutine  is  maintained  for  use  with  the  Naval 
Postgraduate  School,  Monterey,  California,  IBM  360  computer. 

It  was  programmed  by  J.  R.  Ward  in  October,  1963;  was  revised 
in  November,  1971;  and  was  redocumented  in  October,  1973. 
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